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Abstract 

We study the homogeneous Dirichlet problem for the doubly nonlinear equation u t = 
A p u m , where p > 1, m > posed in a bounded domain in M. with homogeneous bound- 
ary conditions and with non-negative and integrable data. In this paper we consider the 
degenerate case m(p — 1) > 1 and the quasilinear case m(p — 1) = 1. We establish the 
large-time behaviour by proving the uniform convergence to a unique asymptotic profile 
and we also give rates for this convergence. 
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We are interested in describing the behaviour of non-negative solutions of the homogenous 
Dirichlet problem for the doubly nonlinear equation (DNLE) for large times. To be precise, 
we consider the following initial and boundary value problem 

{ut(t, x) = A p u m (t, x) for t > and x G tt, 
u(0,x) = uq(x) for x G ft, 

u(t, x) = for t > and x G dfl. 

for to > 0, p > 1. The problem is posed in a bounded domain Cl G with initial data no > 0, 
uq G L l (Q) so that the solution u(x,t) > too. The p-Laplacian operator is well-known to 
be defined as A p w := div(| Vw\ p ~ 2 Vw). We study the large time asymptotic behaviour of 
solutions to Problem |1.1| in the "degenerate case" m(p — 1) > 1, also known as slow diffusion 
case, and in the "quasilinear case" m(p — 1) = 1. 

Let us first make some comments concerning the range m(p — 1) > 1. When p = 2 we 
recover the porous medium equation (PME) u% = Au m with m > 1 while, when m = 1, we 
recover the degenerate p-Laplacian equation (PLE) ut = A p u with p > 2, both well known 
equations in the literature. Notice that in this paper we only require m(p — 1) > 1, that 
also includes cases where either m < 1 or p < 2. The PLE and the PME, as prototypes for 
degenerate diffusion, enjoy many common properties, such as finite speed of propagation and 
the existence of some special (self-similar) solutions, which play an important role in describing 
the asymptotic behaviour for general initial data. In this paper we complete the panorama 
by analyzing in detail the large-time properties of the degenerate DNLE, which combines the 
difficulties of both equations and offers some new challenges. 

The quasilinear case m(p — 1) = 1 is also interesting to study since it inherits some common 
features of the Heat Equation, Ut = Au (which can be recovered when m = 1 and p = 2): 
this equation is invariant under scalar multiplication, and it is known that a general solution 
converges after rescaling to one of the (stationary) solutions of the eigenvalue problem for 
the p— Laplacian operator. However, when (m,p) ^ (1,2) differences appear at the level of 
regularity and qualitative behaviour. While solutions of the HE are C°° smooth, solutions of 
the DNLE have limited regularity due to the degenerate (singular) parabolic character of the 
equations at the level u = (see Fig. [I]). 

The remaining "fast diffusion case" m(p — 1) < 1 has quite different properties and deserves 
a separate study. Indeed, we deal in this case with singular diffusions, and new phenomena 
appear such as extinction in finite time, or lack of uniqueness of the asymptotic profile. All 
this gives a different flavor to the analysis of the asymptotic behaviour. 
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Figure 1: Ranges of parameters m and p 



As references for the previous theory for the DNLE we mention [TO] for the degenerate and 
quasilinear cases and |18j for the singular case. We mention also that the asymptotic behaviour 
of the Cauchy problem on R N has been studied in [TJ. Many of our results are new even in the 
p-Laplacian case m = 1, p > 2. We also remark that most of the techniques needed to prove 
existence, uniqueness and other basic properties of the parabolic DNLE flow, can be taken 
from the books [2D G2] for the PME, and [9] for the PLE. We also refer to [3] and [20] for 
a complete asymptotic analysis of the Dirichlet problem on bounded domains, for the PME 
when m > 1. 

Presentation of the main results. The purpose of this work is to analyze completely 
the asymptotic behaviour of the DNLE on Euclidean bounded domains. For convenience we 
assume that the boundary dQ is C 2,a smooth. Since the cases m(p — 1) > 1 and m(p — 1) = 1 
involve different techniques, we will present them separately. 

Ia. The degenerate case m(p — 1) > 1. This work generalizes the asymptotic analysis 
carried out in the above mentioned papers j3j [20] . The outline of the theory is similar but the 
double nonlinearity asks for a number of interesting techniques. Throughout the study we will 
fix the notation /j, = \/{m(p — 1) — 1) > 0, since this quantity will appear frequently. 

The asymptotic behaviour is better understood via the well-known method of rescaling and 
time transformation; let us introduce 

(1.2) v(t,x) = t^u(t,x), t = e T . 

In this way, Problem |1.1| is transformed into 

v t (t, x) = ApU m (r, x) + fj, v(r, x), for r G R and xeSl, 

( L3 ) S v(t, x) = 0, for r G R and x G dQ, 

v(0,x) = vq for 

In Section [2] we prove Theorem |2.1[ which shows uniform convergence of the rescaled solution 
v(r,x) to its unique asymptotic profile f(x), as r — > +oo. The stationary profile / can be 
characterized as the positive solution to the corresponding stationary problem 

A„/ m + nf = 0, in n, / = on 30. 
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The result of this Theorem is not surprising, but it does not appear explicitly in literature and 
it is needed to prove the next results. The techniques used in this step follow the work (2U] for 
the PME. 

In Section [3] we prove sharp rates of convergence of v(-, r) — > f as r — > 00; this represents the 
first important result of this paper. 

Theorem 1.1. (Weighted rate of convergence) Let Q C M> be a bounded domain of class 
C 2 ' a , a > 0. Let u(t, •) be the weak solution to Problem 1.1 corresponding to a nonnegative 
initial datum uq G Then for every to > fixed there exists C > such that the following 

inequality holds 

(1.4) |(1 + tyu(t,x) - f(x)\ < Cf(x)(l + t)~ 1 for allt>t and x eft, 
where C depends only on p, m, N, uo, O and to. 

In other words 

(1.5) u(t, x) = r»f(x) + 0(f(x) t-^- 1 ). 



Remarks, (i) Sharpness. Estimate (1.4) is sharp, since it is satisfied with equality when we 
consider the separate variable solution, 

(1-6) U(t,x;s) = (s + trf(x), 

with parameter s > 0. 

(ii) Convergence in relative error. Let U(t,x) := U(t,x;0) be the separate variable solution 



(1.6) and let 



v(t, x) = t m( -p- 1 '>- 1 u(t,x) , 



be the rescaled solution given in (2.14). We can rewrite (3.1) in the following form: 



Corrolary 1.1. (Convergence in relative error) Under the assumptions of Theorem \l.l\ 
if u denotes the solution of Problem \l.l\ the 



len 



(1.7) 



u(t, 



U(t, 



L°°(n) 



v(t,-) 



/(•) 



< Ct 



-1 



The main idea will be to compare the general solution u of Problem 1.1 with functions 
belonging to special families, more exactly self-similar solutions of the DNLE. We will try to 
follow the strategy used in the papers [3] and [20] for the case m > 1 of the PME, and solve 
the problems caused by the nonlinearity of the p-Laplacian operator. 

Ib. Selfsimilar study. In the process of proving the above results we became interested 
in the existence and properties of self-similar solutions of the DNLE, that is, functions of the 
form 

U(t,x) = {t + s)- a h{r), r = \x\(t + s)~ p , 
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where a, (5 are positive parameters and h a real valued function satisfying a certain ODE. As a 
by-product we give a formal characterization of such solutions. Selfsimilar solutions are often 
used as a way of indicating the behavior of a general solution of the DNLE. 

II. The quasilinear case. In Section[5]we study the asymptotic behaviour of solutions of the 
DNLE when m(p — 1) = 1. Our study uses the preliminary work [16] and requires a delicate 
barrier technique inspired from the work [6] on fast diffusion stabilization. To be precise, we 
consider the rescaling 

(1.8) v(t,x) = e Xlt S{x), 

where Ai is the first eigenvalue of the p-Laplacian operator A p , [21 [15]. Notice that here there 
is no time transformation, but only rescaling. 

This problem was previously studied by Manfredi and Vespri in [16J, where the authors 
obtained the convergence, along time subsequences, of v(t,x) to a possible asymptotic profile. 
At the same time they proved that the set of asymptotic profiles is included in the set of 
solutions of the corresponding elliptic problem 

(1.9) - A P V = XiV^ 1 inU, V = on dtl, 

and moreover V > in Q. It is known that when A = Ai then the set of solutions is a linear 
set, i.e., they have the form {cV\ : c > 0}, where V\ is a particular normalized solution (a 
normalized p-eigenfunction) , cf . [21 [15] . 

In this work, we complete the asymptotic analysis by proving uniform convergence of the 
rescaled solution v(x,t) to an unique asymptotic profile; this happens for all times t — > oo and 
we also prove a relative error version for this convergence. Our main result in the quasilinear 
case is the following. 



Theorem 1.2. Consider m(p — 1) = 1. Let Q C M> be a bounded connected domain of class 
C 2 ' a , a > 0. Let u(t,-) be a weak solution to the Dirichlet Problem |l.l| corresponding to the 



nonnegative initial datum uq G L 1 (17). Let v(t,x) 
constant c > such that 



u(x,t). Then there exists a unique 



(1.10) 



lim 

t— >oo 



u(t, 



U(t,-) 



1 



L°°(n) 



lim 

t— >oo 



v(t,-) 



o. 



L°°(f2) 



where U(x, t) = e" Al *S(x) and S m = c*/. 



In order to clarify the result, let us point out that the result states that there is a unique 
asymptotic profile of the form S = V l / m , where V is one of the positive solutions of (1.9). 
In other words, there is a unique c = c(uq) > such that V = cV\ =: cf m . Though the 
asymptotic constant c depends on the data, there is no explicit or semi-explicit formula to 
compute it. This is a typical occurrence issue of nonlinear evolution problems, see a discussion 
of the issue in |14] when studying the Barenblatt equation for elastoplastic filtration, a quite 
different model of nonlinear heat flow. In order to prove Theorem |1.2[ the methods used in 
the degenerate case do not work anymore and therefore we apply a different method, a barrier 
argument, based on the one used in [6] to prove convergence in relative error for the fast 
diffusion equation. 
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2 Asymptotic behaviour for m(p — 1) > 1 



2.1 Preliminaries 

In order to present the asymptotic behaviour, we first need to introduce some preliminary 
results concerning the smoothing effects of the DNLE. 

Notations. 

Qt = (0, T) x O, Q = (0, oo) x n, d(x) = dist(x, dn). 

The notion of weak solution is defined in the standard sense, we refer to [9]. In addition, it 
is known by standard semigroup theory that there exists a unique non-negative weak solution 
u of Problem |1.1| with good regularity properties and satisfies the Maximum Principle. 

Now, we illustrate specific properties concerning the smoothing effects of the DNLE, proper- 
ties that will be needed in our proofs (we refer for example to |21j). In what follows we assume 
that O C M> is a C 2,a domain, for some a G (0, 1). Let u(t, •) be a weak solution to Problem 
1.1 corresponding to the nonnegative initial datum uq G L 1 (f]). 



1. Benilan-Crandall type estimates 

1. If m(p — 1) > 1, then 

(2.1) u t > -pd^u 
in the sense of distributions. 

2. If m(p — 1) < 1, then 

(2.2) ut < lit~ x u 
in the sense of distributions. 

Also, in the case m(p — 1) > 1, the weak solution u of Problem |1.1| verifies the following 
estimate 

(2.3) \\ut(t + s,x)\\i <m(* + s) -1 ||«(»)Hi- 

This inequality is a property of viewing the solution u(t) of Problem |1.1| with initial data no 
as a semigroup u(t) = T t (uo). As an immediate consequence we observe that 

(2.4) lk(M)||i < M* _1 ||«o||l 



In addition, using estimate (2.3) with t/2 instead of t and s one can obtain that 

(2.5) K(*,aO||i <^ _1 ll«(*/2)Hi- 
2. Smoothing effects 

In the case m(p — 1) > 1, the solution u of Problem 1 1 . 1 1 satisfies the following 

(2.6) \\u(t, -)\\ L r m < Ct~», t G (0, +oo), r > 1. 
Combining the previous estimates we obtain the absolute bound 

(2.7) KV)||i,oo (n) < Cr^t G (0,+oo). 



As a consequence, we can improve the estimate ( |2.3[ ) by using (2.7) and we get 
(2.8) \u t {t,x)\ x < Cnr l {t/2)~^ < Ct- 1 -^. 
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2.2 Asymptotic behaviour: Uniform convergence to the asymptotic profile 



Now we are ready to state the first important result of the paper. 



Theorem 2.1. Let £1 C M be a bounded domain of class C > a ,a > 0. Let u(t,-) be a weak 
solution to Problem 1.1 corresponding to the nonnegative initial datum uq G L 1 ($7). Then there 
exists a unique self-similar solution of Problem \1.1\ 



such that 
(2.9) 



U(t,x) = t~^f(x), t G (0,+oo),x G n, 



lim t^\u(t,x)-U(t,x)\= lim \t^u(t, x) - f(x)\ = 0, 

t— >+oo t— >+oo 



unless u is trivial, u = 0. The convergence is uniform in space and monotone non- decreasing in 
time. Moreover, the asymptotic profile f is the unique non-negative solution of the stationary 
problem: 



(2.10) 



\f m {x) + ii fix) = o, igsi, f(x) = o, x g on. 



Proof. 1. The main tools are the a-priori estimates (2.7) and (2.1) 

(2.11) u(t,x) <Ct~ M , 
and 

(2.12) u*(*,x) > -C[ir x u. 

These estimates make sense also in the classical way since the weak solution u of Problem |1.1| 
is in fact locally Holder continuous. 

2. Rescaled orbit and equation. We use the rescaling 

(2.13) u(t,x) = t-*v(T,x), t = e T . 

As we have seen, Problem |1.1| is mapped into the rescaled problem: 



(2.14) 



v t (t, x) = A p v m (T, x) + fiv (r, x), for r G M and 

u(0, a;) = fo(x) = u(x, 1), for x G 

f (r, x) = 0, for r G M and x G 5f2. 



For this problem, we take zero Dirichlet boundary data in the sense that v(t,x) G W c \' p (n). 
The initial data are taken non-negative and integrable in The possibility of delaying the time 
origin and the regularity theory allow us to assume that v(x, 0) is bounded, even continuous. 

3. Convergence. The advantage of the new variable is seen when we translate the estimate 
information in terms of v. From the a-priori estimates (2.11 ) and (2.12 ) we get better estimates 
for the function v. 



(2.15) 

and 

(2.16) 



< v < C, 



v T > 0. 
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We conclude from this that for every x £ £1 there exists the limit 

lim v(t, x) = f(x) 



and this convergence is monotone non-decreasing. Also, from (2.15), the function f(x) is 
nontrivial and bounded. Moreover, by (Beppo Levi's) Monotone Convergence Theorem we 
have 

v(t, •) — > f strong in L 1 (S1). 

Since there is an L°°(£l)— bound the convergence takes place in all L q (£l), 1 < q < p. This 
function / is called an asymptotic profile and we will prove that it is the solution of the 



stationary elliptic problem associated to the rescaled problem (2.14) and it is unique. 



4. Energy estimates. We consider the next Lyapunov functional, called entropy, 



EM = ElvM] :=- [ \Vv m (T,x)\ p dx —fi f 

P Jn m + 1 J n 



H I v m+1 (T,x)dx. 



We compute the entropy dissipation (Fisher information) 

d f 

—E(r) = -m / v m ~ 1 Mx)v 2 T {T,x)dx < 0, 
dT Jn 

which means that E{t) is a non-increasing function. Then we obtain that 

1 



(2.17) 
(2.18) 



m 



E{t) = - / \Vv m (T,x)\ p dx 

pJn m + 1 

< E(0) = ^J \Vv^(x)\ p dx - ^PT^ J v™ +l {x)dx. 



-H / v m+1 (T,x)dx 
In 

m+l I 



From this energy estimate (2.17) together with a-priori estimate (2.15) we get that \Vv m (t 



is uniformly bounded in time in every L 9 (fi)-norm, 1 < q < p, thus weakly convergent up to 
subsequences. Let us denote by M the uniform bound for the L p (f2)-norm: 



/ \Vv m (T,x)\ p dx < M, Vt G M. 
Jn 



(2.19) 
We deduce that 

(2.20) d Xi v{r, ■) — v d x J{-) when r — > oo in L q (£l), for every 1 < q < p. 

We can also get a uniform bound for ||w r ||£i(Q), using the Benilan-Crandall type estimate 
(2.8) as follows: 



\\V T (T,X) - flv(T,x)\\ L1{n) = \\A p V m (T,x)\\ L l {n) = || Ap (t^u{t,x)) m [| L l( n) 

= t 1+ »\\A p u m (t,x)\\ L i m = t 1+ ^||n t (t,x)|| L1(n) < d, 
and now we obtain a uniform bound in time. Thus, by now we have: 
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and then 



(2.21) \\v T \\ L i {n) <C 2 . 

Remark: in the PLE case these estimates can be obtained also in the L P (Q) because of the 
contractivity property of the p-Laplacian operator in LP(Q). 

5. Convergence in measure of gradients. The weak convergence of Vv m to Vf m can be 

improved, and in fact we will prove a stronger result, the convergence in measure (with respect 
to the Lebesgue measure C). The idea comes from [5J. We refer also to [7j where the authors 
prove this result for the fast p-Laplacian equation. A strong argument in our proof are the 
following inequalities for vectors in W 1 . If 2 < p then 

(2.22) <\a\ p - 2 a-\b\ p - 2 b,a-b>>~ /1 \a-b\ p , Va, b G TSL N , 
where 71 = c p is a constant depending on p and n. If 1 < p < 2 then 

( 2 - 23 ) < i«r 2 « - i & r% a - b >> 7 2 |a|2 L;~^ |2 . p , 

with optimal constant 72 = c p = min{l, 2(p — 1)}. For a proof of these inequalities we refer to 

USE!- 

We prove now that Vi> m (r, •) converges in measure to V/ m (-) when r — > 00. We remark that 
it is sufficient to prove that (Vu m (r, -))r>o is Cauchy in measure. Thus we have to prove that 
for every e > there exists a T > and A > such that 

£({x£n-. \Vv m (n,x)-Vv m (T 2 ,x)\ > A}) < e, Vn,T 2 > T. 

Let e > and S be the set whose measure we want to estimate 

5 := {x £ n : \Vv m (n,x) - Vt; m (r 2 , x)\ > A}. 

Then 

Sc{xeQ: |Vw m (ri,2;)| > A} U {1 £ SI : |Vw m (r 2 , x)\ > A} 

U{x G n : \Vv m (T U x)\ < A, \Vv m (T 1 ,x)\<A, \Vv m (n, x) - V« m (r 2 , x)| > A} 
= S 1 \JS 2 . 

Since |Vv m (r, -)| is uniformly bounded in L P {£1) then C{S\) < e for 17, t 2 sufficiently large. 
Now, in order to estimate C(S 2 ), the idea is to use algebraic inequalities (2.22) and (2.23) for 
the vectors Vv m (Ti) and Vw m (r 2 ). 

• If p > 2 then 

S 2 c{xefi: \Vv m (n,x)\ < A, \Vv m {n,x)\ < A, 

(\Vv m (T 1 )\ p - 2 Vv m (T 1 ) - \Vv m (T 2 )\P- 2 Vv m (n)) ■ (Vv m (n) - Vv m (T 2 )) > 7l A p := ft}. 

• If 1 < p < 2 then 

S 2 c{x £ £1 : iVt/^x)! < A, \Vv m (T U x)\ < A, 

QVv^h^Vv^n) - \Vv m {r 2 )r 2 Vv m (r l )) ■ {Vv m (r l ) - Vv m (r 2 )) > l<i^r v =: ft}. 
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We remark that in the particular case of the degenerate PLE we only have to consider the case 
p > 2. 

Now, for /3 = /3i if p > 2, respectively j3 = /3 2 if 1 < P < 2, we obtain that 
1 



C(S 2 ) 



s 2 P 
1 

1 
1 

1 
1 

1 



(\Vv m (T 1 )\ p - 2 Vv m {T 1 ) - \Vv m (T 2 )\ p - 2 Vv m (t 1 )) • (Vv m (n) - Vv m (r 2 ))dx 
[V.{\Vv m {T l )\ p - 2 Vv m (r 1 ) - \Vv m (T 2 )\ p ~ 2 Vv m (T 1 ))] [v m (n) - v m (r 2 )] dx 
(A p v m (n) - A p v m (T 2 ))(v m (n) - v m {r 2 ))dx 
(«r(n) + M«(Ti) - Vr(r 2 ) - Iiv(t 2 )) {v m (n) - v rn (T 2 ))dx, 

where we used integration by parts. We recall that v is positive and uniformly bounded in 
time by (2.15) and the norm ||w T ||i2m) is also uniformly bounded in time by (2.21). Thus 

where C is a constant positive number and it follows that 

C{S 2 ) < e, 

for j3 big enough. 

Thus, we proved that the sequence (Vf m (r, -))r>o is Cauchy in me asur e, thus it converges in 
measure to a function W : ft — > Mr. It is well a known fact (Lemma 6.1 in the Appendix) that 
if a sequence is uniformly bounded in L P {Q) and converges in measure, it converges strongly 
in any L q (Q), for any 1 < q < p. It follows that 

Vv m (r, •) -»• W(-) in (L q (n)) N when r -»• oo, for every 1 < q < p. 



Since we already know the weak convergence (2.20) we get that w = V/ m a.e. in fl and we 
can conclude that 

Vv m (T, •) -> V/ m (-) in measure 

and 

(2.24) Vt) m (r, •) V/ m (-) in (L«(J2))^, for every 1 < g < p. 
Thus, we get that, up to subsequences, 

(2.25) Vv m (T, •) -»• V/ m (-) a.e. in 



6. The limit is a stationary solution. Multiply equation (2.14) by any test function 
4>(x) £ C%°(Q) and integrate in space, x E fl, and time between n and r 2 = ri +7b, for a fixed 
T > 0. We get that 

r r~i~2 r rri r 

v(r 2 )4>dx — / v(Ti)4>dx = I / A p v m (j)dxdt + fi I / vipdxdt = 
Jn Jn Jn Jt x Jn 



n Jn 



\Vv 



Vv m V(j)dxdt + fj, 



Tl Jf7 



vcbdxdt. 
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Let t\ — > oo. Then also T2 — > 00 and we get that 



v{T2)4>dx — I v{T\)(j)dx — > 
n Jn 

and then, by Lebesgue's Dominated Convergence Theorem, it follows that 

T2 



fi / v<j)dxdt — > Tofi / f(f>dx. 
J 71 Jn Jn 

Now, we need to obtain the convergence of the integrals involving gradients: 

(2.26) ( 2 [ \Vv m \ p ~ 2 Vv m V<j)dxdt T I |V f m \P- 2 V f m Vcpdx 

Jn Jn Jn 

when n — > 00 in order to obtain that 

= -T / |V/ m | p - 2 V/ m V^x + / uT / 
and dividing by To and integrating by parts we get 

= J A p f m 4)dx + n J f4>dx, 



which proves that / is a weak solution of the stationary problem (2.10) 

-A p f m (x) = fif(x), x^n. 



In order to justify assertion (2.26), we remark that, after a change of variables, the left term 
can be written as 

[ T ° [ \Vv m (T + T 1 ,x)\ p - 2 Vv m (T + T 1 ,x)V<pdxdt. 
Jo Jn 

Thus, it is enough to prove that 

f ° I \Vv m (T + n,x)\ p - 2 Vv m (T + n,x)V(f>dxdt^T [ \V f m \ p - 2 V f m V(pdx when n -> 00. 
Jo Jn Jn 



The idea will be to use Lemma |6.2| from the Appendix in the following context. Let Q± = 
[0,To) x 17 (finite measure space), H = ~S, N (Hilbert space) and let 

A:H^H, A(Z) = \Z\ P ~ 2 Z. 

Then A is single valued, monotone and R(A + I) = H and therefore, according to the theory 
of monotone operators, A is maximal monotone. We consider the sequences 

Z n (t, x) = Vv rn (T + n,x):n 1 ^ H, 

W n {T,x) = A(Z n {r,x)) = \Vv m (T + n,x)\ p - 2 Vv m (T + n,x) : fli -> H. 



11 



The hypothesis of Lemma |6.2| are satisfied as follows. From (2.25) 

Z n (T,x) -> Z(t,x) = Vf m (x) a.e. on 



Let q = Then W n (t, x) is uniformly bounded in L q {VL\] H), by the energy estimate (2.19), 
since ^ 

f ° [ ||Vw m (r + n, x)\ p ~ 2 Vv m (T + n, x)\ q dxdt < MT , Vn > 0, 
Jo Jn 

and thus it converges weakly (up to subsequences) to a function W in L q (£l\; H) when n —> oo. 
But since Q± is bounded, then weak convergence in L g (Qi;H) implies the weak convergence 
in L 1 (f2i;ff) and thus 

W n (T,x) W(r,:c) weakly in L x (Oi;iJ). 



Now, according to Lemma (6.2), we obtain that 

W(r,x) = A(Z(r,x)) = |V/ m (x)r 2 V/ m (x). 
Thus, the weak limit of W n {j,x) is unique and we get that 

\Vv m (T + n,x)\ p - 2 Vv m (T + n,x) -» |V/ m (x)| p ~ 2 V/ m (x) weakly in L«(fi i; ff). 
By taking smooth and compactly supported in Q we obtain the desired convergence 

f ° [ \Vv m (T + n,x)\ p - 2 Vv m (T + n,x)V$dxdt^ [ ° [ |V / m | p_2 V f m V <\>dxdt 
Jo Jn Jo Jn 

= T [ \Vf m \ p - 2 Vf m V4>dx. 
Jn 

7. Uniqueness of the stationary solution. Let us prove that the nonnegative and nontrivial 
stationary solution is unique. If we have two stationary solutions of (2.10), f\ and fa, we can 
construct solutions of the DNLE of the form 

U 1 {t,x) = r tM f 1 {x), U 2 (t,x) = (t + s)-»f 2 (x), 

for some s > 0. U2 has initial data U2(x,0) = s~ ,J 'f2(x). Formally, U\(x, 0) has infinite values 
and then by the Comparison Principle we conclude that U2(t,x) < Ui(t,x). The technical 
details of the proof are as follows: by the L 1 —dependence theorem of weak solutions of Problem 
11.11 we know that 

' (C/ 2 (t,x)-f/i(t,x)) + ^, 



is decreasing in time. The proof of this fact is standard: we perform the difference of the 
equations of U\ and U2, multiplying by h(w), where w = U™ — U™ and h is a C 1 (]R) function 
such that < h < 1, h(s) = for s < and h'{s) > for s > 0, and then integrate on 0. 
Notice that 

< /J-o := inf h'(w(x)) < 00. 
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Because of the nonlinearity of the p— laplacian operator, we will use again the algebraic in- 
equalities p32j) and (pT23). 



(U 2 (t,x) - Ux(t,x)) t h{w)dx = / (A p U^(t,x) - A p f7f (t, x))h{w)dx 

n 



(|VC/n p - 2 VC/r - \VU^\ p ~ 2 VU^) Vh(w)dx 

= - [ {\VU^\ p - 2 VUr - |VC/ 2 m | p - 2 VC/ 2 m ) ti(w)V(V? - U?)d2 
Jn 

< -vo [ {\vu^\ p - 2 vu^ - \vu^\ p - 2 vu^) v(t/r - U^)dx. 

Jn 



Now, if p > 2 we obtain that 

(U 2 (t, x) - Ui(t, x)) t h(w)dx < -/i 7l 

i 

and if 1 < p < 2 the estimate will be 

(U 2 (t,x) - U\(t,x))t h(w)dx < -/io72 



|VC/f - VU^\ p dx < 0, 



WW 



VC/ 2 m | 2 



|VC/f p-f + |vc/ 2 m | 2 -p 

Letting /i converge to the sign function signQ we get that 

d 



dx < 0. 



dt 



(U 2 (t,x) - Ui(t,x))+dx < 0. 



Now, the integral goes to as t — > because U\{t,x) goes pointwise to infinity as t — > and 
then U 2 (t,x) > Ui(t,x) for t large enough. We conclude that (U 2 (t,x) — Ui(t,x))+ = a.e. in 
x for every t > 0. Thus U 2 (t,x) < Ui(t,x) a.e. in x for every t > 0. Using the form of Ui and 
U 2 , we get 



/ 2 (x) < 



t + s 



Letting s — > we get f 2 (x) < fi(x). The converse inequality is similar. 



8. Better convergence. We have established the result (2.9) in the sense of con- 
vergence. The passage to uniform convergence depends on having better regularity for the 
solutions , i.e. on a compactness argument. As we mentioned before, (6.1), uniformly bounded 
solutions of the DNLE are C a continuous in space and time with uniform Holder exponent 
and coefficients. 

Consider now the second type of rescaling that we may call fixed-rate rescaling 
(2.27) u\(t,x) = X^u(Xt,x). 

For every A > the function u\ is still a solution of the DNLE to which the a-priori estimate 
(2.11) applies. Hence, in a set of the form (1,2) x Q this family is equi-continuous and by 
Ascoli-Arzela Theorem it converges along a subsequence Xj — > oo. Now, observe that 

u\(l,x) = v(logX,x) 

to conclude that v(log Xj, x) converges uniformly. Since the limit is fixed, /, the whole family 
v(t, x) converges as r — > oo and (2.9) is proved. □ 
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2.3 Some remarks on the asymptotic profile 



1. Existence. The proof of Theorem 2.1 also guarantees the existence of a solution of the 
stationary problem (2.10) by obtaining it as the limit of v(t, •) when t goes to oo. As we have 
previously established, this solution is called the asymptotic profile of parabolic problem ( 1.1 ). 

Furthermore, we recall a second proof of existence, based on an entropy method, which can 
be applied also for the general case of solutions with changing sign. 



The elliptic problem (2.10) can be written in terms of the function w = f m (the notation 
makes sense since / > in f2 by Maximum Principle) as 



(2.28) 



A p w(x) + [iiv™ (x) = 0, x G £2, w(x) = 0, x G c%7. 



The typical approach to solving equation (2.28) for the experts in elliptic equations is to view 
the solution w as a critical point of the functional 



(2.29) 



J(g) 



1 

PJsi 



Vw\ p dx 



m 



m + 1 



m + l 

w m dx. 



The proof is classical and we resume it following the ideas from |20j . It can be showed that 
J is well defined in Wq ,p {£V) since m(p — 1) > 1, J is bounded from below via Poincare's 
Inequality, and also the infimum is negative. Moreover, along any minimizing sequence there 
is convergence in Wg' p (Jl) and the infimum is taken, hence J has a minimum. Also, 

J(g) > J(w), Vg G W^ p (fl), 



where w is the solution of (2.28 ) and it follows that w is the point where J attains its minimum. 

2. Uniqueness. We already proved the uniqueness of the asymptotic profile in point 7 of the 
proof of Theorem |2.1| 



3. Regularity. We know that w is a bounded solution of equation 

Apw(x) + (j,w™(x) = 0. 

By known regularity results for degenerate elliptic equations, we get that w G C^(SV) for some 
6(0,1]. 

4. Behaviour near the boundary. Concerning the behaviour of /, the following estimates 
were proved in |16j : 

(2.30) |V/| < C d{x) 1/m ~\ Vx G fl, 
and 

(2.31) dd(x)^ m < f(x) < C 2 d(x) 1/m , Vx G 0. 

that is w(x) = f m (x) has a linear growth near the boundary. 
Also, w satisfies a Boundary Principle 

dw 



(2.32) 



dv 



(x) = Vw(x) ■ v{x) < 0, for x G d£l, 



where v{x) denotes the outward unit normal vector to dil at the point x. Notice that, in the 
case of the PLE, the boundary principle is satisfied by /. 
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3 Rate of Convergence for m(p — 1) > 1. Proof of Theorem 



1.1 



In this section we will give the proof of Theorem We previously showed in Theorem 2.1 
that for any initial datum no G L 1 (0), the corresponding rescaled solution t^u{t,x) converges 
to an unique asymptotic profile / uniformly in space and monotone nondecreasing in time. 
The goal of this section is to provide sharp convergence rates, namely to prove that 

(3.1) |(1 + t)»u(t,x) - f{x)\ < Cf(x)(l + ty 1 for ah t > t , x G 0, 

where / is the solution of the elliptic problem ( 2.10[ ) 

A p / m + nf = 0, x G n, f = on dQ. 



The proof of this result is based on the techniques introduced by Aronson and Peletier for 
the PME in [3j. Although a similar proof can be adapted to the case of the DNLE with some 
lengthly arguments, in this work we will give a simpler proof based on the results of Section 
[2] Let us explain the strategy of the proof. 



1. Improved upper bound. In Theorem 3.1 we prove that there exists a constant si > 
depending only on p, m, N, uq and O such that 

o < u(t,x) < (si + fr A 7(x), Vsefi, t>i. 



2. Positivity. In Proposition 3.2 we prove that even if uq has compact support there exists 
T" > depending only on p, m, N, uq and £1 such that 

u(t,x) > 0, Vx G SI, t > T'. 



3. Sharp lower bound. In Theorem 3.2 we prove that there exist T" > and sq > 
depending only on p, m, N, uq and S7 such that 

u(t, x) > (s + t)~' x f(x), VxGfi, t > T". 



Then, the estimate (3.1) follows as a consequence of the upper and lower bounds together 



with the boundedness of the asymptotic profile < / < C. 



3.1 Reduction 

Let f2 C ~K N be a a bounded domain of class C 2,a , a > 0. Let u(t, •) be a weak solution to 
the Problem |l.l| corresponding to the nonnegative initial datum no £ L 1 (17). Previously, in 
Theorem |2.1 1 we proved that the rescaled solution t^u(t,x) is monotone increasing in time and 
convergent to the function f(x), thus bounded from above by f(x): 

(3.2) u{t, x) < r»f(x), Vt G (0, oo), x G H. 

In virtue of this result we can assume that the data satisfy the following conditions, denoted 
as Hypothesis (H): 

(HI) $1 C WL N be a a bounded domain of class C 2,Q , a > 0. 
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(H2) no is a nonnegative function denned on Q such that no € L 1 ($7), no = on dQ and there 
exists si > such that 

(3.3) uq(x) < s\f(x), Vx € 17. 

Since the DNLE is invariant under time displacement then (H2) is satisfied by starting with 
initial data u(to, •), where to > 0. We assume henceforth that such a displacement in time has 
been done. 

3.2 Improved upper bound for u 



In the following theorem, we will improve the upper bound (3.2) of n that we have previously 



proved in Theorem 2.1 



Theorem 3.1. (Quantitative upper bound) Assume that Q and uq satisfy the hypothesis 
(H), and let u be the corresponding solution of the Problem \l.l\ Then there exists a constant 
s\ > such that 

(3.4) u(t, x) < (si + t) -/ 7(z) Vt>0, xeO. 

where s% depends on p, m, N, uq and fL 



Proof. The proof of (3.4) relies on the Comparison Principle for the DNLE. Consider as com- 



parison function the separate variable solution U of the DNLE given by 

U(t,x) := U{t,x;s 1 ) = {s 1 + t)- fl f{x), 



where s\ > is a constant given in the Hypothesis (H2) which satisfies the inequality (3.3) 

uq(x) < U(0,x) = sif(x), Vx G £1. 
Therefore by comparison, it follows that 

u(t,x) <U(t,x), Vt > 0, xefL 

□ 

3.3 Positivity of u 

In this subsection we prove the positivity of n in Q, under the Hypothesis (H). This will be 



done in two steps. First, in Proposition 3.1, we will prove the positivity of n in a domain 



fli s C fl and then we complete the result by proving positivity of n up to the boundary in 



Proposition 3.2. In this direction, we will make use of the properties of the distance to the 



boundary function d(x) = d(x, d£l) stated in Subsection 6.3 of the Appendix. In terms of d(x) 
we define the following sets 



Mir = {x E O : d(x) > r}, Q r = O \ = {i £ Q : d(x) < r}. 
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Proposition 3.1. (Inner positivity) Assume that £1 and uq satisfy (H) and let u denote 
the weak solution of Problem \l.l\ and f denote the solution of problem (2.10). Let < 25 < £o 
fixed, where £o is defined in Lemma\6.3\ Then there exist e > and T\ G [0, +oo) such that 



u(Ti,x) > e for all x G 
where e and T\ depend only on m,p, N, Q and uq. 



Proof. The main tools are the uniform convergence (Theorem 2.1) of the rescaled solution 

v(t, x) = t^u(t, x), t = e T , 



defined in (2.13) to the asymptotic profile / and the properties of / given by (2.31). More 



Let e = 



exactly, there exists C\ , C 2 > such that 

(3.5) Cidm(a?) < fix) < C 2 d™(x), Vx G SI. 

> 0. Then there exists Tq > such that 

\\v{t,x) - f(x)\\ L *c(p,) < e , Vr > T . 
Since the convergence of v(t, x) to f(x) is monotone nondecreasing in r we derive that 

v(t,x) > f(x) — eo, Vx E f2, Vr > Tq, 



and then, by using (3.5), for a; G Qj^ we obtain the lower bound 

v(t, x) > Ci<im (x) — eo > Ci^m — eo = eo- 
In terms of u(t, x) these estimates rewrites as 

u(t,x) > e t~^, Mx £ n I<s , Vi > e T °. 

Let 



(3.6) 
Then 



e = e rr M = 6 e- T ^. 



u(Ti,x) > e, Vx G ^7,5. 



□ 



Proposition 3.2. (Positivity up to the boundary) Assume that f2 and uq satisfy the 



hypothesis (H) and let u denote the weak solution of the problem 1.1 and f denote the solution 



of the problem (2.10). Consider T\ > given by Proposition 3.1, Then there exists T2 > 

u(Tx +T 2 ,x) > 0, \fxen, 



and T 2 depends only on m,p, N, f2 and uq. 
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Proof. We consider < 25 < £o as in Proposition 3.1 First, by (6.4), we observe that 



^28 C [J B 2S (y). 



Then, since we have already proved the positivity inside the domain in Proposition 3.1 it is 
enough to demonstrate that there exists T 2 > such that 

u(T 1 +T 2 ,x) > 0, VxG B 2S (y), Vy G dSlj^g. 

Let e given by (3.6). Let y G d£li j2 s and consider the Barenblatt solution U as in Section [4] 
such that 

supp U{x — y, 0; a, s) = Bg(y) and max U{x — y, 0; a, s) = e, 

that is 



(3.7) 



0(m(p-l)-l) / CX m(p-l)-l 

' 8 = 5 L 



Now, assume the time T 2 = t when supp U(x — y, t; a, s) riches the boundary of f2, that is when 

supp U(x-y, t; a, s) = B 2S {y) 

which implies 



(3.8) 



s. 



We want to apply now the Parabolic Comparison Principle. To this aim we need to compare 
u(T\ + t,x) and U{x — y,t;a,s) when (t, x) belongs to the parabolic boundary {0} x B 2 $ U 
[0, T 2 ] x dB 2 $. Firstly, for t = and x G B 2 $, we have 

.. (T s > (e>U(x-y,0;a,r), x G B s (y); 
u K n,x)_^ = u{x _y j0 . ajT)j x£ B 25 (y)\B 5 (y). 

Secondly, when t G [0, T 2 ] and x G dB 2 $ we have 

u{T\ + t,x) > 

Therefore we obtain that 

u(Ti +T 2 ,x) =U(x-y,Ti +T 2 ;a, s), Vx G B 25 {y) 



e >U(x — y,0;a, r), x G dB 2 $ \ dQ; 
= U(x - y,0;a,r), x£dB 2S C\dn. 



We notice from (3.7) and (3.8) that T 2 does not depend on the point y, but only on the data. 

□ 
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3.4 Sharp lower bound for u 



In the following theorem we will derive a lower bound for u, similar to the upper bound we 



have previously proved in Theorem 3.1 



Theorem 3.2. (Quantitative lower bound) Assume that f2 and uo satisfy the hypothesis 
(H) and let u be the weak solution of the problem and f be the solution of the problem 



(2.10). Then there exist two positive constants sq > and T4 > such that 
(3.9) u(t,x) > (so + tyvfix), Vx G H, Vt G [T 4 ,+oo), 

where so and T4 depend only on m,p,N,U and uo. 

Before we start the proof of this theorem, we will establish the following preliminary results. 
Let 

V(x - y,t; M, s) = (t + s)- a [g^(\x - y\(t + s)^)} + 
be the so called intermediate selfsimilar solutions defined in Section [H 



Proposition 3.3. Under the assumptions of Theorem 3.2, for every y G dQj : 25, where < 
25 < £0 is fixed as in Proposition \3.1\ we can choose constants M and s such that there exists 
a time T3 > ; for which the next inequality holds: 

(3.10) u(Ti + t,x)>V(x-y,t;M,s), \fx G U, Vt G [0, T 3 ], 



where T\ is the time we have obtained in Proposition 3.1 and T3 > depends only on m,p, N, Q, 
and uq and it is independent of y. 



Proof. The same ideas as in Proposition 3.2 apply since the functions V have a similar behaviour 



as the Barenblatt functions. Consider the selfsimilar subsolution V such that 
supp V(x — y, 0; M, s) = B$(y) and max V(x — y, 0; M, s) = e, 

that is 

M = es a , 8 = S Pa(M) = /M m(P p 1)_1 a(l) and s = ( e -(™(p-i)-i). 



o(l) 

The exact values of M and s are not important; what matters is that they depend only on the 
data and, in particular, they are independent of y. 

Now, consider the time T3 = t when supp V(x — y, t; M, s) riches the boundary of J7, that is 
when 

(3.11) supp V(x-y,t;M,s) = B 2S (y) 
from where we deduce the explicit value for T3 

/ 25 \ 1/(3 

(3.12) T 3 = I 2 *) -s. 
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Then the Parabolic Comparison P rinci ple can be applied to u and V on the parabolic domain 
[Tx,Tx + T 3 ] x H as in Proposition 



3.2 



and we obtain that 



u(Ti+t) > V(x-y,t;M,s), Vx G O, Vt G [0,T 3 ]. 



We notice from ( |4.9| ) and (3.12) that T3 does not depend on the point y, but only on the 
data. □ 



We define 
(3.13) 



T1+T3 



where T\ and T4 are given by Proposition |3.1[ respectively Proposition 3.3 



Proposition 3.4. (Boundary behaviour) Under the hypothesis of Theorem 3.2, let T4 as 



in (3.13) and 5 as in Proposition 3.1, Then there exists a constant uj > 



u m (T 4 , x) > ujd(x) for igSIj. 
such that uj depends only on m,p, N, $7 and uq. 



Proof. Let y G ■ Then by Lemma (6.3), there exists a unique z(y) G dfl such that 

S = d(y) = \y- z(y)\. Let 

V(x - y,T 3 ; AT, a) = (T 3 + s)- a [<7-(|z - y\(T 3 + 
be the self-similar subsolution obtained in Proposition |3 . 3| where M and s are given by formulas 



(4.9). Let [0, a) the largest interval starting from where g > and by (3.11) it follows that 
a = a(M) = 2S(T 3 + Then, in view of the continuity of g', there exist ko < k\ < such 

that 

ko < g'(v) <ki, V77 G [5(T 3 + s)"' 3 , 2<5(T 3 + s)~% 

and it follows that 



(3.14) 



9(V) > \ki\(a- V ), Vr? G [S(T 3 + s)~ p , 25(T 3 + 



Thus it follows from (3.10) and (3.14) that for every iGfion the segment between y and z(y) 
such that 5 < \x — y\ < 25, that is d(x) < 5, u can be bounded from bellow as follows 



(3.15) 
where 



n m (T 4 , x) > V m ((x - y, T 3 ; M, s) = (T 3 + s)' 

> M(t 3 + s )— m (2<5(r 3 + s )-^ - 

= |fc 1 |(T 3 + S )^ m ^(25-|x-y|) 
= ujd(x), 



*g(\x-y\(T 3 + s 
■-y\(T 3 + s)- 



UJ 



\h\(T 3 + s) 



G M + is a constant which depends on the data, but not on y 
and x. Observe that (3.15) holds for arbitrary y G d£lj 25 and for all x on the inward directed 



normal through y provided that d(x) < 5. As we remarked in Lemma (6.3), the normal map 

1 r G [0,£o)- Therefore, it follows 

u m (T4,x) > ud(x) for x G Qs- 



H r is a homeomorphism for all r G [0,£o)- Therefore, it follows from (3.15) that 
(3.16) 



□ 
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Proof of Theorem 3.2 First we will prove there exists k 2 > such that 
(3.17) u(T 4j a;) > k 2 f(x), VxGO. 



By Proposition 3.3 and Theorem |2.31| tt satisfies 

u m (T 4 ,x) > ud(x) > ujC^f m {x), for igSI,. 



Moreover, by Proposition |3.1| and the boundedness of the profile < / < C in Q we obtain 
that 

u(Ti,x) > e > — fix), VxeU ItS . 

Then, since T\ < T4, inequality (3.17) is satisfied with k 2 = min{e/C, w 1//m C2}. 

Finally, let U(t,x) = (so + be the separable solution of the DNLE with initial data 

s o M /> where so is defined by the relation Sq M = k 2 . Then 

u(T 4 ,x) > U(0,x), Vx G 0,, 



and (3.9) follows by Comparison Principle. 
End of the proof of Theorem S.2\ 



4 Study of self-similar solutions for the DNLE 

In this section we will illustrate a short description of the self-similar solutions of the DNLE, 
focusing on the properties useful for the proofs in the paper. A complete analysis of these 
solutions is beyond the purpose of our paper and it can be the subject of a future work. 

For a complete characterization of self similar solutions of the PLE in the case p > 2 we refer 
to [5]. Likewise, for the relation between self-similar solutions of the PLE and those of the 
PME we make reference to |12| . 

Self-similar solutions of the DNLE are functions of the form 

U(t,x) = {t + s)~ a h(r), r = \x\(t + s)~P, 
where s > is a constant, a and f3 are positive parameters related by 

(4.1) (m(p- 1) - l)a + pp = 1. 

The profile g := h m : [0, 00) — >• K is a function satisfying the differential equation 

(4.2) ag^(ri) + /3r (g^j (r) + ^l|^( r )|^V(r) + (p - l)\g' {r)\^ 2 g" (r) = 0, r > 0, 
also written in the equivalent form 

(4.3) a/ l (r) + /3r/i / (r) + -^ rT (r iV - 1 | 9 , (r)| p -V(O)) / = 0, r > 0. 

Self-similar solutions are (possibly signed) solutions of the DNLE in the whole space. When 
the support of the positive part of such a function U is included in f2 then IA + is a sub-solution of 
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the DNLE equation. For this reason, self-similar solutions are useful to indicate the behaviour 
of a general solution u of the DNLE. 

We consider the initial conditions 
(4.4) h(0) = M m , h'(0) = 0. 



The existence of a positive solution of ODE (4.2) with initial conditions (4.4) on an interval 
[0, a), where a € (0, oo], can be proved using fixed point methods when p < 2 and using phase 
plane methods when p > 2 (we refer to [3] when p > 2 where the author discusses the case of 
the PLE). As far as we know, fixed point methods do not work when p > 2. 

If one multiplies the ODE by r^" 1 and integrate between and r, where r £ [0, a), it follows 
that 

(4.5) \g'(r)r 2 g'(r) = -f3rg^(r) - f ' s N - l g^(s)ds, 

r Jo 

or, equivalently, 

(4.6) |^ (r)r V(r) + ^=^ s N( g iy {s)ds= ® £( r) . 

J 

We will make a formal study of self-similar solutions of the DNLE by considering the following 

cases: (3 = 0, a — (3N = 0, a — j3N > and a — (3N < 0. We define the numbers 

1 a OLE 1 1 

a B '■= — 7 7\ -i , / /An i p B := — = - — — — , a := 



m{p - 1) - 1 + (p/N) ' (m(p-l)-l)JV + p' m(p- 1)-1" 




(a) Case /3 = 0, a = qq (b) Case a = j3N 

Figure 2: Self-similar solutions of the DNLE 

I. Case p = 0, a = oo. Separate variables solutions 

They have the form 

u(t,x) = (t+ s rf(x), 



where s > is a constant and / is the solution of the stationary problem (2.10). Notice 
that the functions belonging to this family are self-similar solutions according to the previous 
definitions when / is a radial function. These functions are very useful since they indicate the 
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asymptotic behaviour of u, the general solution of the DNLE and thus we will use them for 
comparison. A second aspect is that they do not propagate and thus we have to consider also 
different types of self-similar solutions. 

II. Case q = (3N. Barenblatt solutions 

For more details we refer to [211. A Barenblatt solution (also called source type solution) 
exists for the DNLE in the "good range" 

P 

m(p - 1) + — > 1 

that includes of course m > 1 (the degenerate PME) and p > 2 (the degenerate PLE). When, 
moreover, m(p — 1) > 1 (our case), Barenblatt solutions have the form: 



U(x,t;a,s) = c(t + s)~ a (aA - \x(t + s) -/3 |A) 



m(p — 1) — 1 



where s > is a positive parameter and 
(4.7) 

1 cub ( m(p - 1) - 1 / a 

a = a B = — 7 7T — : ■ , , P = Pb- " 



m(p — 1) — 1 



m(p - 1) - 1 + (p/N) ' N ' \ p \N 

When s = 0, this function has a Dirac delta as initial trace 

lim U(x,t) = MS (x). 
t— >o 

The remaining parameter a > is free and can be uniquely determined in terms of the initial 
mass / Udx = Mq. 



Barenblatt solutions are compactly supported and they propagate with finite speed. We will 
use them as a lower bound in order to prove the positivity of u inside fi. Since they have a 
flat landing contact (zero derivative at the boundary of their support) we can not obtain a 
quantitative lower bound for u up to the boundary of £1. Their advantage is that they have an 
explicit formula which is very advantageous for computations. 

III. Case a > f3N. Intermediate self-similar solutions 

In this case 

a > as, < (3 < fis- 

This family of self-similar solutions, which we denote by V, inherits some useful properties of 
the Barenblatt solutions and the separate variables solutions: V has a compact support that 
propagates and g = h m has a transversal cross through the r axis. This is explained as follows. 
Consider [0,a) the largest interval starting from where h > 0. Then, from (4.5) we obtain 
that g'(r) < for r G [0, a) and thus 

\g'(r)r 2 g'(r)<-Prg^(r). 

Furthermore, this implies that 

i i i 
—g (r) > /3p-i rp-i g m (p-i) . 



23 



Integrating from to r with g(0) = M m we obtain that 

m(p — 1) 

/ m(p-l)-l 77J,(l3 — 1) — 1 1 P \ m(p-l)-l 

5(r) < M p- 1 — p P -i r p-i ] , 

\ mp ) 

for all r € (0, a). From this upper bound we derive an estimate for the point a where h(a) = 

p-i 

m(p-l)-l / mr) \ p 1 

a<M p — ^ pp. 

\m(p — 1) — 1 J 

The important fact is that a is finite, thus V has a transversal cross through the r axis at the 
point r = a with 



(4.8) 9'{a) = -(°^- P s N -V /m (s)ds 

\ a Jo 



i 

' 1 =: k < 0. 



We will use the form 

V(x,t;M, S ) = (t + s)- a [h(r;M)} + , 
and therefore the following characterization of the support 

supp V(x, t; M, s) = {(x, t) : \x\ < a(t + , t > 0}. 

We denote by 

a = a(M), k = k(M), h(-) = h(-;M) 

in order to emphasize their correspondence to the Cauchy problem with initial conditions 
h(0) = M, h'(0) = 0. We remark that 

m(p— 1) — 1 m(p— 1) — 1 

(4.9) h(a(M);M) = 0, h(r;M) = Mh(M p r; 1) and a(M) = M p a(l). 

Subsequently we will consider the self-similar solutions described above only on [0, a), the 
largest interval starting from where they are positive; for complete definition, on [a, oo) they 
are assigned zero values. This way, they are sub-solutions of the DNLE. 

Remark. In the present work we do not study the behavior of these functions when then g 
takes negative values. Depending on the values of a and /?, g can behave differently, as we can 
see in Figure 3. 



(a) Case a> /3N (b) Case a> /3N (c) Case a < f3N 

Figure 3: Self-similar solutions of the DNLE 
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IV. Case a < /3N 

The self-similar function corresponding to the profile g in this case does not have a compact 
support, hence this class is not useful for our estimates. For completeness we will provide a 
formal characterization of these functions. 

Recall that in this case 

< a < o>b, (3 > /3b. 

Using basic computations as in the previous case, one can easily prove that g is a positive 
decreasing function converging to when r goes to oo. Moreover, for every tq > there exists 
C = C(r ) > such that 

am 

(4.10) g{r) >Cr P , Vr > r . 

Asymptotic decay. We point out that g behaves as r — > oo like G a {r) = ar -7 , 7 = am/ (3. 
We continue with a formal proof. 

We will consider the following series expansions of g: 

g(r) = a\r _1 + ... 

as r goes to 00, where 7 > is an exponent to be determined and is representative for 
lower order terms. Then 

ah(r) + /3rh'(r) = af (a - ft— \ r~™ + 

V mJ 

pvU {r N - l \9'(r)r 2 g\r)))' = a?~V'~ 1 (7(P - 1) +p- N)^ 1 ^ 1 ^ 1 + .... 

For a comparison of the first terms we notice that 

_( 7 + l)(p_l)_l<_2. 

m 



Thus the leading asymptotic term in the expansion of the ODE formula (4.3) is 

af (a - [3— ) r~™. 
V mJ 

Moreover, we notice that the coefficient a — (3j/m = 0, from where we deduce that exponent 
of the leading asymptotic term is 

(4.H) 7= T - 

At this time we have no information about a\. The coefficient of the remaining term is 

al'Y-H^p-lJ+p-N), 

whose sign depends on the values of f3. 
Let 

1 1 

Pb ■= 7 7 T~\ 1 m\ »r ' a B '■ 



(m(p — 1) — 1 + p/N) N ' m(p- 1) - 1+p/N' 

N-p m(p-l) 1-pPi 
71 := r> Pi '■= l~i 7\ 1 1 /An at = m{p-l)(3 B , «l :^ 



p-1' (m(p — 1) — 1 + p/N) N ^ ' r ' m(p-l)-l' 

The sign of coefficient ^{p — \) + p — N can be now obtained depending on the values of (3. 
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1. Case j(p -l)+p-N>Oo 



2. Case j(p - 1) + p - N < 



7 > 71, (3 < Pi, a> ax\ 
or p > N. 

P G (fii,l/p), a<ai, p<N; 
impossible, p > N. 



3. Case 7(j>— — A^ = 44> 7 = 71, (3 = fii, a = ot\. Then <7(r) = r 7 is the solution 

of the ODE (Ob. 



We can observe that in the first two cases we can not deduce the decay of g and we need to 
perform a second approximation. 

5 The quasilinear case m(p — 1) = 1 



In this section we consider Problem ( 1.1 ) for m > 0, p > 1, posed in a bounded domain Q, £ 
with smooth boundary of class C 2,a and initial data uq > 0, no € L 1 (0). We study the large- 
time asymptotic behaviour of solutions of the problem in the quasilinear case m(p — 1) = 1. 

As usual, the problem is better understood via the method of rescaling. We consider 
(5.1) v(t, x) = e xt u(t, x), t € [0, 00), x E CI, 

where A is a real parameter whose choice we will justify in the next subsection. Then v is a 
solution of the rescaled problem 

{v t (t, x) = A p v m (t, x) + Xv(t, x) for t > and x e £1, 
v(0,x) = Uo(x) for x£(l, 

v(t, x) = for t > and x G d£l. 

The study of the asymptotic behaviour in the present case m(p — 1) = 1 differs considerably 
from the case m(p — 1) > 1 previously studied for several reasons. Firstly, the proof in the case 
m{p— 1) > 1 is based on monotonicity: the rescaled solution t^u{t, x) /* v{t, x). This argument 
cannot be applied in the present case. Besides, we have no universal a-priori estimates similar 
to the degenerate case. 

5.1 The associated stationary problem 

Consider the stationary problem associated to Problem |5.2[ 

(5.3) A p f m + Xf = in Q, f(x) = for x G dQ, 
that can be rewritten in terms of z = f m as 

(5.4) A p z + X\z\ p ~ 2 z = in ft, z = on dfl. 

We want to obtain solutions z(x) > in Q. We call eigenvalues the A-s for which there 
exists a nontrivial solution of Problem |5.4[ which is known as the eigenvalue problem for the 
p— Laplacian. The following result was proved in [2| 115]. 
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Theorem 5.1. (Simplicity and isolation of the first eigenvalue of Problem 5A) The 

first eigenvalue X\ of Problem |5.4| is simple and isolated. Moreover, Ai is the unique positive 
eigenvalue of Problem |5.4| having a nonnegative eigenf unction. 



Moreover, the first eigenvalue Ai of Problem |5.4| can be characterized as 



that is Ai = C~ p where C is the best constant of the embedding W 1,p (fi) into L P (Q). 



(o.5) A i :-m\-\ / p G < p (fl) 



5.2 Preliminary estimates for the evolution problem 

We state two results obtained by Manfredi and Vespri (Theorems 1.4 and 1.4 from |16j). 

Theorem 5.2. Consider m(p — 1) = 1. Lei C 1^ 6e an open, bounded domain. Then 
there exists a unique solution of the problem \l.l\ corresponding to an initial datum uq G L 1 (J1). 
Moreover, for all t > 1, there exists a constant c(t) such that 

(5.6) \u(t,x)\ < jie~ Xlt c(t)f(x), VxGfi, 



where f is a solution of the problem 5.3 such that f G C°(f2), / P_1 V/ G L p (£l) and 71 is a 
positive constant depending only on the data N, p, m, the L 1 norm of uq and the C l,a norm 

ofdn. 



Theorem 5.3. Consider the hypothesis of the previous theorem and assume moreover that 
uo > and not identically zero. Then, for every t > 1, there exist the constants c(t), c(t) G M+ 
such that the following estimate holds 

(5.7) e- Xlt c(t)f(x) < u(t,x) < e- Xlt c{t)f(x), Vx G fl 
Moreover we also have 

(5.8) 7i(t)e _Al *d(a;) p - 1 < u(t, x) < j 2 (t)e- Xlt d(x) p - 1 , Vx G O, t > 1 
and 

(5.9) |Vn(t,x)| < j 3 (t) e - Xlt d(x) p - 2 , Vx€Sl,t>l. 



Remarks 



Reduction. The estimates given by the previous theorems are true for every t > to, 
where to > is fixed. Since the doubly nonlinear equation is invariant under a time 
displacement, we can assume that the previous estimates are valid for every t > 0, 
otherwise we can start with initial data u(to). We assume therefore such a displacement 
in time has been done. 

We can fix / in any way, up to a multiplicative constant. Therefore we fixe / a nonneg- 



ative solution of the problem 5.3 
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Inspired by the ideas from |14j we will obtain more information about the constants c(t) and 
c(t). Let us define 

(5.10) c(t) = inf{c : v(t, x) < cf(x)}, c(t) = sup{c : v(t, x) > cf(x)}. 



According to the Theorem |5.3| the previous definition makes sense: c(t) < oo and c(t) > 
for every t > 0. Thus, we can take c(t) and c(t) to be the best constant such that estimate 



(5.7) holds. It is a simple consequence of the Maximum principle that c(t) and c(t) are two 



decreasing, respectively increasing functions of t. Therefore the following limits exist: 

(5.11) Coo = lim c(t), c(t) \c~oc, 

t— >oo 

(5.12) c 00 =limc{t), dt)^^. 



In addition, the constants c(t) and c(t) are uniformly bounded 
(5.13) C < c(t) < < Coo < c(i) < d, Vt>0. 

We can sum up what we have proved so far in the following lemma. 



Lemma 5.1. Let v be a solution of the rescaled problem (5.1). Then there exist the positive 

constants Co, Ci, C2 > such that 

(5.14) Qoo f{x) < v(t, x) < Coo/(x), Vx G n, t > 0, 

(5.15) C d(x) < v(t, x) < Cid(x), Vx£n,t>0, 
and 

(5.16) \Vv(t,x)\ < C 2 d(x) p ~ 2 , \/xen,t>0, 



where Co,Ci,C2 > depend on Q, Ai and f is the positive solution of problem 5^3_ we have 
taken. 

As a consequence we obtain the uniform convergence, up to subsequences, of v(t, •) to a 
stationary profile. 

Theorem 5.4. (Uniform convergence to an asymptotic profile up to sequences) 

Consider m(p — 1) = 1. Let $7 C ~R N be a bounded domain of class C 2,a , a > 0. Let u(t, ) 
be a weak solution to Problem \l.l\ corresponding to the nonnegative initial datum uq G L 1 {Q). 
Then for any given T > there exists a sequence r n — > 00 such that 



(5.17) \e Xl ^ +t ^u(T n + t,s)-c*f(x)\ ->0, r n 



00, 



uniformly for x G and < t < T ; where f is the positive solution of problem \5.3\ we have 
taken and c* is a positive constant. 
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Proof. I. Energy estimates 

la. We consider the following energy functional 

E(t)=E[v(t)] :=- f \Vv m (t,x)\ p dx — Ai / v m+1 (t, x)dx. 

P Jn m + 1 



n 



We compute the energy dissipation 

d 



E(t) = I(t) = m I v m - x vldx > 0, 
at ./n 



which means that E{t) is a non-increasing function and 

-t 2 



E{t l )-E{t 2 )= r i{t)dt. 



As well, we deduce that the integral 

rt 

I(t)dt 



is convergent as t — > oo and E(t) has a limit as t — > oo. 
Since v(t, x) is bounded in SI uniformly for t > we obtain that 

(5.18) [ \Vv m (t,x)\ p dx < M, Vt>0, 

Jn 

in other words |Vt> m (t, -)| is uniformly bounded in L p (£l) for t > 0. 



Ib. As a consequence of (5.18) one can prove via Holder's Inequality the following technical 
result 

(5.19) [ (A p v m (t l ,x))v m (t 2 ,x)dx < M, Vti,t 2 >0. 

Jn 

II. Convergence. We define 

(5.20) v T (t,x) = v(t + r,x), t,r>0,xeO,. 



Then v T is still a solution of Problem 5.2 with initial data {v(0,;c) = v(t, x). 



We fix T > 0. The family {v T ) is relatively compact in 

X = L°°([0,T] x O) 

thus it converges along subsequences 

v Tn (t,x) — )• S(t,x) uniformly in (t, x) € [0, T] x 
From the a-priori estimates we deduce the boundedness of v 

C < v(t,x) < Ci, Vt>0, xGQ, 
and since S is the limit of v Tn (t,x), then S also satisfies the same lower and upper bounds. 
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In what follows we fix such a subsequence (r n ) and the corresponding limit S(t,x) 



III. Convergence in measure of gradients 

Similar to the case m{p — 1) > 1 one can prove the convergence in measure of the sequence 
(Vt>™(-, •))„, where in the present case 

v Tn :[0,T]xn^ [0,oo). 

More exactly, we can prove by similar methods that the sequence (Vt>™) n >o is Cauchy in 
measure, thus it converges in measure to a function W : [0, T] x Q, — > M. . It is a well known 



fact (Lemma 6.1 in the Appendix) that if a sequence is uniformly bounded in L p and converges 
in measure, then it converges strongly in any L q , for any 1 < q < p. It follows that 



Vu™ — > W strongly in (L q ([0,T] x fi)) when r n — > oo, for every 1 < q < p 
Thus, we get that, up to subsequences, 



(5.21) 

and we conclude that 



V<( v )->W(v) a.e. in [0,T] x fi, 
W(t,x) = VS m (t,x). 



IV. The limit is a solution of the stationary problem 



Multiply equation (5.2) by any test function 4>{x) E C£°(0) and integrate in space, x E f2, 



and time between r n and r n + T. We get that 

'■r n +T 



(5.22) / {v{T n +T)-v{T n ))4>dx 



rTn+T 

Vv m \ p ~ 2 Vv m V<pdxdt + X 1 / 

Q J r n 



vcj)dxdt. 



(i) The left hand side term of (5.22) is uniformly bounded independently of T: 



(5.23) 



v{r n + T)(pdx - I v{T n )(j)dx 
a Jn 



<2Ci|n|||0|| L oo ( n). 



Furthermore, if <j) is supported in a compact K C where < c\ < v < C2 and < s < T 
then 



(5.24) 



/ (v{T n + s) - v(r n )) 4>dx 
Jn 



< 



\vt(t) 4>\dxdt < 



1/2 



< r':/- 1 2 f / U f {{v^l 2 ) t fdxdt\ =CT 1 ' 2 



I(t)dt 



1/2 



Since the double integral J x I(t) is finite, the integral on the right hand side goes to zero as 
r n — > oo. On the other hand, 



/ v(r n + s)4>dx- / v (r n )(j)dx = / v Tn (s, x)4>dx - / v Tn (0, 
Jn Jn Jn Jn 



x max 



S(s,x)(f)dx— / S(0,x)4>dx, r n — > oo. 
n Jn 
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Thefore, we showed that 

r (S(s,x) -5(0, x)) <j>dx = 0, 



n 



for any < s < T and any test function <p with supp (j) = K CC 0. Since if is arbitrary, the 
result holds for all (j) G C£°(0) which implies that 5 is independent of time on [0,T]: 

(5.25) S(s) = S(0), VsG[0,T]. 

(ii) For the right hand side, we continue as follows. Let r n — > oo. Since Vw™ — >• W = VS m 



a.e. in [0, T] x O then using Lemma 6.2 of the Appendix we get that 

T n +T p pT p p 

/ \Vv m \ p - 2 Vv m V(j)dxdt = / / \Vv™J- 2 Vv™V(i>dxdt^T |VST~ 2 V5 m V<^x- 

io is! Jn 

The last integral 

[■Tn+T P j-T I- 

j j v(j)dxdt = / / v Tn (t,x)(f>(x)dxdt — >■ T / S(t,x)4>(x)dx, r n — > oo. 
Jt„ Jn io Jn Jn 

Therefore 

(5.26) - / |VS m | p - 2 V5 m V^x + Ai / S(pdxdt = 0, Vte[0,T] 

Jn Jn 

and thus S is a weak solution of the stationary problem |5.3| According to the Theorem 5.1 

S(t,z) = 5(0, x) = c*/(x), Vt G [0,T], x G 0. 
For simplicity, we denote := c*f(x). 

□ 

Remarks 

• V = S m , where S := c*f, is a positive solution of the eigenvalue problem for the p— Laplacian: 

(5.27) A P V + XiVP- 1 = in fi, F = on 311 



• As a consequence of (5.11) and (5.17) the constants c^, Coo and c* satisfy 

Qoo — ^* — ^oo ■ 

• The previous proposition does not guarantee the uniqueness of a stationary limit S. There- 
fore the rescaled solution v(t, x) may oscillate between the bounds c 00 f{x) and Coof(x) by 
converging on subsequences to asymptotic profiles of the form c*/ with < c* < Cqq. This 
kind of behavior has to be considered in the case of some parabolic evolution equations, for 
example in the case of signed solutions of the porous medium equation ut = Au m , m > 1. The 
set of possible asymptotic profiles is obtained as the w-limit of the solution and it is contained 
in the set of classical solutions of the associated stationary (elliptic) problem (we refer to the 
survey [20]). 

Next, we will prove that an oscillating behaviour is not possible. More, exactly, we show that 

Qoo = C* = Cqo, 

which guarantees the existence of a unique asymptotic profile S = c*f and therefore the 
uniform convergence of the rescaled solution v(t, x) to S for all times, uniformly in S7. 

To this aim, we will study the behaviour of the quotient v/S up to the boundary. 
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5.3 The relative error function and its equation 



Assumptions (A). In what follows we make the assumptions: C IR^ is a bounded domain 
of class C 2,a , a > 0, u(t, •) denotes the weak solution to Problem in the case m(jp — 1) = 1 
corresponding to the nonnegative initial datum uq G We fix T > 0, the corresponding 



sequence r n — > oo and the constant c* 6 [c^jCm] obtained in Theorem 5.4 for which the 
rescaled solution v(t,x) = e Xlt u(t,x) converges 

[| e Ai(rn+*) u ( Tn +t ,.)- c*f{-)\\ L oo m -)■ 0, r n oo 



uniformly for i £ [0, T], where / is the positive solution of the problem 5.3 we have taken. We 
denote S := c*/ and we call it possible asymptotic profile. 

Starting from this partial convergence result, we will obtain a much stronger convergence as 
t — > oo: we show the uniqueness of the asymptotic profile and the convergence in relative error 
of v(t, •) to S(-) up to the boundary as a consequence of the next proposition and estimates 



(5.14). 



Proposition 5.1. (Behaviour up to the boundary) Under the assumptions (A) there exists 
a unique constant c* > depending on uo and Q, and for given e > there exists t(e) > such 
that 

v m (t x) 

(5.28) ~ 6 < s ^ x) J — 1 < e, \fxeCl, Vt>t(e). 



Moreover, S = c*/ is i/ie solution of the problem 5.3 announced in Theorem 5.4 



Motivated by the techniques used by Bonforte, Grillo and Vazquez in [6] we will use the so 
called relative error function and the method of barriers. 

To this aim, we introduce the Relative Error Function(REF) 
v m (t r) 

(5.29) ( t, x) = _l^_Z_l, v m = S m^ + 1)= y^ + 1)j and V = S m, 

Notations. We define 



= {x £ £1 : d(x) > 5}, Qg = O, \ Q It s = {x G f2 : d(x) < 5}. 



where in what follows 5 > is considered to be a small positive parameter (see Subsection 6.3 
of the Appendix for properties of the distance to the boundary function). 

Properties of the REF 

• The parabolic equation of the REF. Using the equations satisfied by v and V and relation 
m(p — 1) = 1 we obtain that 

(5.30) (p - 1)(1 + 0)^ = V^ p -^A p ((0 + l)V) + Ax(0 + l)*" 1 . 



• is uniformly bounded in (t,x) for t > 0. This can be derived from the estimates (5.15) on 
v and 5, which is a stationary solution: 



, , i 1 — C2.m < < C3 j?Ti , 



32 



In any interior region fljs C SI, the REF function satisfies 



1 + 



V 



> in r2/ 5< 5 for any t > 0. 



Regularity of solutions of the parabolic equation (5.30). Since (f> is also bounded in the interior 



of ft, we conclude that the parabolic equation (5.30) is neither degenerate nor singular in the 
interior of U. Also the solution <f> of such a parabolic equation is Holder continuous in any 
inner region C fl since both v and S are Holder continuous and positive in the interior of 

n. 

Convergence of the REF in an interior region of ft. Under the running assumptions, 



we know by Theorem 5.4 that 



sup|?;(T„ + t) 
n 



as n 



oo, 



uniformly for t £ [0, T] for a fixed T > and a corresponding sequence (r n ) n , but this is not 
sufficient to prove the convergence of the quotient v m /S m to 1 in the whole Q, since at the 
boundary there is the problem caused by the fact that both v and S are and therefore the 



parabolic equation (5.30) may degenerate at the boundary. However such a problem is avoided 



in any interior region where both v and S are strictly positive. 
We can sum up the results we proved so far in the following lemma. 



Lemma 5.2. (Inner convergence) Letv the solution of the reseated problem (5.2) and S the 
solution of stationary problem (|5.3l) corresponding to a given T > and a sequence r n — > oo 



as in Theorem 5.J±. Let (f> be the associated relative error function defined by (5.29). Then 

-+0, 



as t„ 



ijn + t, ')\\L°°(n I>s ) = SUP |0(T„ + t, •] 

oo uniformly in x S Qj $ and < t < T, for any given 5 > 0. 



5.4 Construction of the upper barrier and consequences 



In order to finish the proof of Theorem 1.2 we have to prove the uniform convergence of <fi up 



to the boundary. This will be realized using a barrier argument based on the ideas from [6]. 
Throughout the paper we will use the notation £o f° r the critical value which implies good 



properties of the distance to the boundary function in f^ as we explain in Subsection 6.3 of 
the Appendix. 

Let us first point out some connections between distance to the boundary function and the 



solutions of the eigenvalue problem (5.27). 



Lemma 5.3. (Properties of the asymptotic profile V = S m ) Let V be a solution of the 



eigenvalue problem (5.27). Then V satisfies the following estimates: 



1. There exist Co > and C\ > such that 

C^d(x) < V(x) < C?d(x), Vx G ft. 
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2. For every < £1 < £o there exists a constant /3q > such that 

VV(x) • Vd(x) > p > 0, Vx G %. 

3. For every < £i < £o ^ere exists a constant K\ > suc/i i/iai 
(5.31) < liTi < |W| < K 2 , VsG%. 



Proof. Point 1. is a consequence of Lemma 5.1 



The proof of the point 2. is similar to the one given [6] since the function V involved has the 
same properties as its correspondent in the fast diffusion problem. 



Point 3. is a consequence of the estimate from point 1. and estimates (5.16) and (6.3) since 
| W| = ImS^VS] < K 2 d(x) p - 2+ ^ m - 1 ^ m = K^x)^' 1 ^ 1 = K 2 , Vx G fi. 

□ 

Next, we present the construction of the barrier that plays an important role in the estimate 
of REF (j) close to the boundary. We mention that our construction is different from the one 
of [6] where the operator was the usual Laplacian A. In our case, the p-Laplacian operator 
contains also mixed derivatives of second order whose estimate is more technical. 

Lemma 5.4. (Upper barrier) We can choose positive constants A,B,C so that for every 
to > the function 



(5.32) 



$(t, x) = C- BV(x) - A(t - t ), 



is a super-solution to equation ( 5.301) on a parabolic region near the boundary 



= {(t,x) G (t ,oo) : $(t,z) > -1}, 
and moreover C (to, To) x fi^, where £i < £o- 



Proof. We will prove that the function (5.32) is a supersolution for the equation (5.30) on the 



parabolic region if we can find constants A, B and C such that 

(5.33) (p - 1)(1 + $) p - 2 $ 4 > y-( p ~ 1 )A p (($ + 1)V) + Ai(* + l) p -\ 
We will prove that a convenient choice for A, B and C will be of the form 

(5.34) (Ai(C+l) + A(p-l))^ _1 <uB, 
where 



min{l,2 2 - p } • 



2(p~l)K p 



From the beginning we assume that (t,x) G S$ C (to, To) x f2^ 1 where To is such that 

C- BV(x) 



34 



The left hand side term satisfies 

(p - i)(i + $) p " 2 ^ = -(p - + $) p ~ 2 - 

The right hand side term is of the form 

V- {p -^A p ((^ + 1)V) + Ai($ + If' 1 = V-(P-VA p f(V) + Ai($ + l)^ 1 

where 

/(z) = (C + l-Bz-A(t-t ))z. 
The term A p / can be computed as 

\f(V) = \f'(V)r 2 f'(V)A p V + (p- l)\f'(V)r 2 f"(V)\VV\P. 
Properties of the function / 

1. Function / is a concave parabola with zero values at the points z = and z = zo where 

C + l-A(t-t ) 
zo ■= g • 

2. The derivatives are 

f'( z ) = C + l-2Bz-A{t- t ), f"(z) = -2B. 

3. When applied to V, the derivative 

f'(V) = & + l-BV. 

Moreover, sufficiently close to the boundary, f'(V) is positive and bounded. By choosing 

fr or\ C • L 1 Z C+l-A(t-t ) \ 

(5-35) 6 = mm | &> — - = — j, 

we obtain the following bound on : 

< V(x) < C?d(x) < < j 

and then 

o<^</ / (y(x))<fc 1 , 

where 

(5.36) h := /'(0) = C + 1 - A(t - i ), y = /' (f ) = ^(C + 1 - A(t - to)) > 0. 

Since p > 1, we obtain a lower bound for f'(V) p ~ 2 on as follows 

(5.37) f'(V) p ~ 2 > fik p - 2 , := min{l, 2 2 ~ p }. 
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Sufficient conditions for the parameters 

Since V is a solution of the stationary problem (5.27) then — A P V = Ail^ -1 and inequality 



(5.33) can be rewritten as 



-A(p-l)(l + <S>) p ~ 2 > -\ 1 \f'{V)\ p - 2 f'(V)-2B{p-l)V~( p - 1) \f{V)\ p - 2 \VV\ p + \ 1 (<Z> + l) p - 1 . 

The idea is that can have large values close to the boundary, thus it is sufficient to 

find A, B and C such that 



(5.38) 



Ai($ + If" 1 + A(p - 1)(1 + $) p - 2 < 2B(p - l)V-( p -V\VV\ p \f'(V)\ p - 



For £i as in (5.35) and the bounds (5.31), (5.37), the right hand side term of (5.38) satisfies 



the lower bound 

2B{p-l)V^ p ~^\VV\ p \f{V)\ p - 2 > 2B(p — l)(3K p k p ~ 2 (Ct^ 1 )^ =: II. 



-i 



For the left hand side term of (5.38) on we obtain the upper bound 



Ai($ + l) p - 1 +^(p-l)(l + $) p - 2 < Ai(C+l-A(t-t )) p - 1 + A(p-l)(C+l-^(t-t )) p - 2 :=I. 
Thus it is sufficient to take A, B and C such that 

(C + i _ A (t _ to ))P-2 (Al(C + i- A (t- t )) + A(p - 1)) < 2B(p - l)f3K p k p - 2 



According to (5.36) this inequality becomes 



Ai(C+ 1 - A(t - t )) + A(p - 1) < 2(p - l)/?^- 1 -^. 

si 



One can see that a sufficient condition on A, B and C would be 



Ai(C + 1) + A(p - 1) < 2(p - l)/3iff Cf 1 , . 

si 



B 



□ 



We will obtain an upper bound of the REF <j> at a certain time Ti up to the boundary as a 



consequence of comparison of <fi with the barrier function of Lemma 5.4 



Lemma 5.5. Let be the barrier function introduced in Lemma 5.4- given by 

$(t,x) =C- BV(a?)-i4.*, 
Let r n — > oo 6e a sequence along which the REF converges to 1 as stated above. Then for every 



e > we can choose n e > and positive constants A, B, C and 5 as in Lemma 5.4 such that 

(5.39) <f>(t + T n ,x)<$(t,x), VxeSls, Vn>n e , Vt e [0,Ti], 

where 

C - BCf <5 - e 



(5.40) 



ri = Ti(e,<$) 
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Proof. We fixe e > and consider < S < £1 where £1 > is given as in Lemma 5.4 Also, let 
T > and (r n ) that we fixed in Assumptions (A). By the uniform inner convergence stated in 
Lemma 5.2 we know there exists n tt s > such that 

(5.41) \4>(t + T n ,x)\<e for x G Qi,s, t G [0,T], n > n £j(5 . 

Once we will choose 5 > we will obtain n e as above. 
A first condition on the parameters will be that 



(5.42) 



ri(e,*) = ^ <T. 



Now, we consider the barrier function <3? and prove that cj>(t + r n ,x) < <3?(t, x), for a fixed 
n > n e , on the set (0, Ti) x Qg, where T\ = Ti(e, 5). More exactly, inequality ( |5.39 | follows as 
a consequence of the parabolic maximum principle on this set. 

Therefore, we have to check that this comparison is satisfied on the parabolic boundary 
formed by three pieces x fl s U (0, Ti) x dtti tS U (0, Ti) x dfi. 

1. Comparison of <fi with $ ai i/te initial section t = 0. We want to obtain that 

(5.43) </>(r n ,x)<$(0,x) = C-BV(x) 

for all x G £lg- This is possible because of the uniform boundedness of (f) 



Ci 



1 



C 2 , m < < C; 



3,m 



G: 



for all x S n as a consequence of bounds (5.15). Now, we simply choose C sufficiently large 
and A, B to satisfy (5.34). 

2. Comparison of '(f) with <E> on the inner parabolic boundary. This part of the boundary is given 
by the set 

(0,Ti) x dSl IfS = {(t,x) : t G (0,Ti), d(x) = 5}. 

For (i, x) as before, 5>(i, x) is bounded as follows 

C- At- BC^b < $(t, x) < (7 - At - BC^S. 



Let us fixe e > and < 5 < £i where £i > is given as in Lemma (5.4). By (5.41) 

(f)(t,x) < e for x G fij )(5 , i G [0,Ti]. 
Thus one can obtain < $ if 

(5.44) e < C - At - BC^5, Vi G [0, Ti] . 



Since C can not be small, this implies a choice for A and -B compatible with (5.34) from the 
construction of the barrier. This can be realized by choosing 5 > small. Once C and B are 
chosen it is sufficient to take At small. 

3. Comparison of <f) with on the outer lateral boundary. This part of the boundary is given 
by the set [0, Ti] x dQ, where we only know that (ft = v m /S m — 1 is bounded. As in |6] we can 
use an approximation trick using the solutions of problems posed in the domain £l k C Q. 
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We will prove the desired comparison (5.39) for the function Uk and obtain it for u by passing 
to the limit. 

We know that /* u as k — > oo uniformly on the compact set [0, t] x Q, for every t < T\. 
Then 



W 



1 = -1<0 on [0,Ti] x dn k , 



where IA = e Xlt S(x) is a separate variables solutions of the DNLE in Q. Thus by (5.44) we 
have 

<t> k < < C - At = $ on [0, Ti] x dn k . 

Steps 1 and 2 hold also for Uk since < u. Thus, by the parabolic comparison principle we 
obtain that (f>k < ^ on the region fJ fc n fia for i 6 [0, T\]. Passing to the limit when k — > oo we 
obtain cj) < $ on [0, Ti] x O a . 

We obtain in this way an improvement of the upper bound of 4> near the boundary after some 
time delay given by 

C - BC^5 - e 



t<T 1 =T 1 (e,S) 



A 



which is the maximum that (5.44) allows. Notice that the delay time T\{e, 5) does not depend 



on the time r n we fixed at the beginning. 
Therefore, in order to choose the desired parameters we perform the following steps: we 



choose C sufficiently big to have (5.43). Then choose A and B to satisfy (5.34). Finally we 
choose 8 small such that (|5.44|) and (|5.42[) hold , that is t < Ti(e, 8) < T. □ 



Better estimate from above for <fi up to the boundary 



Under the assumptions of Lemma 5.2 and Lemma 5.5 we deduce that for t = r n + Ti(e, 8), 



where Ti(e, 8) is given by (5.40) and n > n e , the REF eft satisfies the upper bound 



(5.45) 



e, d(x) > 8; 

e + BC^8, d\x) < 8. 



Therefore, by fixing e > 0, finding a barrier with constants A, B and C and then taking 
8 < e/(BC™), we obtain the time Ti(e, 8) and the level n e such that for all n > n t we have 



(5.46) 



(r„ + Ti,x)<2e ViGfi. 



This means that v{T\ + r n ) < (1 + e)S. The maximum principle implies now that the 
comparison is valid for all times t > T\ + r ne . This proves that Cqo < , thus they are the 
same. One of the consequences is that c* does not depend on the subsequence, therefore the 
whole family v(t, •) converges to S = c*/ as t — > oo. Moreover, we conclude the uniqueness of 



the profile c*/ as well as the upper approximation stated in Proposition 5.1 



5.5 Construction of lower barriers 

It remains to prove a similar bound for the REF <fi from below. To this aim we define 

v m (t,x) 



jp ■- -(f) = 1 



S m (x) ' 
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We perform a similar approach as in the upper barrier case. 

• The parabolic equation of ip 

(5.47) - (p - i)(i - ^y- 2 A = v-^a p ((i - ^)v) + Ai(i - ipy- 1 . 

• The function tp is uniformly bounded in (t, x) for i > 0. This can be deduced from the 
estimates (5.15) on v and S, which is a stationary solution: 

= C2,m < 1p < C-&,m = 1 - 



1 



Cx 
Co 



Co 



In any interior region Qj g C 0, the function ip satisfies 



1 — if) = — > in 0/ j( 5 for any t > 0. 

Lemma 5.6. (Lower barrier) We can choose positive constants A',B',C so that for every 
to > the function 



(5.48) 



*(t, z) = C + BV(s) - A'(t - t ), 



is a super- solution to equation ( 5.47[ ) on a parabolic region near the boundary 

T 



(t,x) G (t ,oo) : < <Jf(t,x) < 



ant 



moreover C (ioi^b) x Q^ 2 , where £2 < £o- 



Proof. We will prove that the function \J r given by (5.48) is a supersolution for equation (5.30) 
on the parabolic region if we can find constants A', B' and C such that 

(5.49) - (p - 1)(1 - ^) p - 2 ^ t > y-( p - 1 )A p ((l - + Ai(l - ^y- 1 . 

From the beginning we assume that (t,x) G C (io,To) x f^ 2 where To is such that 

(5.50) 



C< + -1/2 C" + 5T(x) 
< Jo — to S - 



A' A' 
The left hand side term of ( 5.49| ) is positive on 

(5.5i) - (p - i)(i - <ay- 2 v t = { P - i)A'{i - my- 2 . 



The right hand side term (5.49) is of the form 



where 
and 



y- (p - 1} A p ((l - V)V) + Ai(l - ^y- 1 = V-^ApgiV) + Ai(l - 

o(z) = (l-C-B2f + A(t-*o))z, 
A P g(y) = \g'(V)r 2 g'(V)A p V + (p - l)|«A^)rVW|W| p - 
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Properties of the function g 

1. Function g is a concave parabola with zero values at the points z = and z = zq where 

l-C' + A'(t-to) 
Z0 '= B> • 

2. The derivatives are 

g '(z) = 1 - C - 2B'z + A(t- t ), g"{z) = -IB'. 

3. When applied to V, g'(V(x)) is positive and bounded sufficiently close the boundary. 
More exactly, we consider 

(5.52) £ 2 = mm <^ £ 



k Cf 4 AB'C" 2 
For this choice we obtain the following bound on 0^ 2 : 

< V(x) < C[ n d{x) < Cf 6 < j 

and then 

0<Y<9'(V(x))<h, 

where 

(5.53) /, ; := ,'(0) 1 C- Ait - to), g' (: ^ ' ( " f ■ VU ^ h) - h 



4 ) 2 2 



Sufficient conditions for the parameters 

Since V is a solution of the stationary problem (5.27) then —A p V = \iVP~ 1 and therefore 
the supersolution inequality (5.49) can be rewritten as 
(5.54) 

(p- l)A'(l -my~ 2 > -\ l \g\V)\ p - 2 g\V)-2B\p-l)V-^-^\g\V)\ p - 2 \VV\ p + \i{l-^)' p - 1 . 

Next, the idea is that when we are sufficiently close to the boundary, y - ^™ 1 ) will be large 
enough and then the right hand side term of (5.54) will be negative. More exactly, on 

-2B'( P - l)V- {p - 1) \g'(V)\ p - 2 \VV\ p + Ai(l - tf)*" 1 < -2B'(p - l)(Cr6) _(p_1) ^2 + Ai, 

where ki is a positive constant given by 

\g' [V)\ p - 2 \VV\ P < max{k p - 2 , (h/2) p - 2 } ■ K p =: k 2 . 

Moreover, if £2 is sufficiently small 

then 

-2B'( P - i)v-( p - 1) \rf(y)\ p - 2 \vv\ p + \ 1 (i-v) p - 1 < 
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and inequality (5.54) holds true. We remark that the above conditions on the distance to the 
boundary £2 can be rewritten as 

C 1-C + A'(t-t ) f2B , (p-l)k 2 ^ 1/(p - ir 



(5.55) 



£2 = min < £ , 



AB'C r { 1 



□ 



Lemma 5.7. Let ^ be the barrier function introduced in Lemma 5.6, given by 

V(t,x) = C + B'V(x)-A't. 

Let r n — > 00 be a sequence along which the REF converges to 1 as stated in Assumptions (A). 
Then for any e > we can choose n e > and positive constants A' , B' , C and 5 as in Lemma 
15.61 such that 

ip(t + T n ,x) < *(i,x), Vx G tt s , Mn > n e , W G [0,T 2 ], 



(5.56) 
where 

(5.57) 



Proof. We fixe e > and consider < 5 < £2 where £2 > is given as in Lemma 5.6 Also, let 
T > and (r n ) that we fixed in Assumptions (A). 

We adapt the proof of Lemma |5.4| by writing the estimates in terms of the function ip = —(f). 
By the uniform inner convergence stated in Lemma 5.2 we know there exists n e > such that 

(5.58) |V>(r„ + t,x)\ < e for x G fi/,<s, * G [0, T], n > n e . 

We impose a first condition on the parameters A', B 1 , C 
(5.59) 



T lM = a + B,C ? S -'<T. 



A' 



Now, we consider the barrier function constructed in Lemma 5.6 and prove that ip(t + 
i~n,%) < tff(t,x), for a fixed n > n e , on the set (0, T2) x Q5, where T2 = T2 je)( 5. More exactly, 
inequality (5.56) follows as a consequence of the parabolic maximum principle on this set. 

Therefore, we have to check that this comparison is satisfied on the parabolic boundary 
formed by three pieces OxfijU (0, T 2 ) x dn IjS U (0, T 2 ) x dVL. 

1. Comparison of tp with at the initial section t = 0. We want to obtain that 

4>{T n , x) < *(0, x) = C + B'V(x) 
for all x G Fig. This is possible because of the uniform boundedness of ip 



C 3 ,m < 1p < C^ m = 1 



for all x G as a consequence of bounds (5.15). Now, we simply choose C > 1. 
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2. Comparison of ip with ^ on the inner parabolic boundary. This part of the boundary is 
given by the set 

(0,T 2 ) x 8Q I)S = {(t,x) : t G (0,T 2 ), x G ft, d{x) = 5}. 
For (i,x) as before, ^(t,x) is bounded as follows 

C - A't + B'C^5 < $(t, x) < C - A't + B'C^S. 



Let us fixe e > and < 5 < £2 where £2 > is given as in Lemma (5.6). By (5.58) 

ip(t + r n ,x) < e for x G i G [0, T2]. 

Thus one can obtain ^ < ^ if 

(5.60) e < C - At + B'CffS, Vt G [0, T 2 ] , 

or equivalently 

«< c ' +fl 7 , -' = ^o. 

3. Comparison of ' tp with ^> on the outer lateral boundary. This part of the boundary is given 
by the set [0, T2] x <9f2, where we only know that ip = 1 — v m /S m is bounded. Here we will 
use an approximation trick using the solutions of problems posed in an extended domain 
ft fc D ft. Like in Lemma 5.5 we prove the desired comparison (5.56) for the approximating 
function Uf~- We know that u k \ u as — > 00 uniformly on the compact set [0,t] x ft, for 
every t < T 2 . Formally 



W 



-00 on [0, T 2 ] x <9ft, 



where IA = e 1 .S'(x) is a separate variables solutions of the DNLE in ft. On the other hand, 
on [0,T 2 ] x <9ft, 

x) = C - A't > C" - A'T 2 = e - BC^S 

and thus the comparison ip^ < holds true on the outer lateral boundary. Finally, steps 1 
and 2 hold also for Uk since Uk>u and thus, by parabolic comparison we obtain that <p\~ < <3> 
on [0, T%] x Qg. Passing to the limit when fc-^oowe obtain ip < ^ on [0, T 2 ] x 

We obtain in this way the an improvement of the lower bound of p near the boundary after 
some time delay given by 



t<T 2 (e,5) 



C + B'CJfS - e 
A? ' 



which is the maximum that (5.60) allows. Notice that the delay time T 2 (e, 5) does not depend 
on the time r n we fixed at the beginning. 

Therefore, in this case of comparison with lower barriers, in order to choose the desired 
parameters we perform the steps: we choose C' > 1 big enough to have the inner comparison 
and then we choose A', B' such that condition (5.60) holds. Finally choose 5 small enough 

□ 



such that condition (5.59) holds. 
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Better estimate from below for <p up to the boundary 



Under the assumptions of Lemma 5.2 and Lemma 5.7 we deduce that for t = r n + ?2(e, 8), 
for a fixed n > n t , where T2(e, 5) is given by (5.57) the REF ip = —(f) satisfies the upper bound 

(5.61) 



e, d(x) > 5; 

e - B'C^S, d(x) < 6. 



Therefore, by fixing e > 0, finding a barrier with constants A' , B' and C and then taking 
5 < e/(-B'C™), we obtain a time T2(e, 5) and a n t such that for all n > n e we have 



(5.62) 



An immediate consequence is the lower estimate of v in terms of the profile v(T2 + r n ) > 
(1 — e)^. The maximum principle implies now that the comparison is valid for all times 
t > T2 + r ne . This proves that > c* , so that they are the same. Therefore we proved that 



This concludes the part of the lower approximation stated in Proposition 5.1 



Once we proved the uniqueness of the asymptotic profile and, Theorem 1.2 follows as a 



consequence of the estimates (5.46) and (5.62). 




6 >J; S 2 S 3 <> ;0.0: 

(a) Behaviour of x) 





*1 


















(2)/ 




d(x) 








2 5 

(b) Behaviour of &(t,x) 



(0,0) 



Figure 4: Idea of the behaviour of the barriers: y-axis: values of <fr(t, x), x-axis: values of d(x) = 
d(x, 90), i.e. the distance to the boundary. Ej: the points where the barrier (a) 4>(i, x) — e,, (b)\&(i, x) = 
ej, i = 1,2, 3. e,: different values of e( decreasing with i=l,2,3) give different barriers 4^, 4^ decreasing 



with e as the arrows (1) and (2) indicates. £1 and 5 as in Lemma 5.4 £2 as in Lemma 5.6 



6 Appendix 

6.1 Two convergence results 

The following lemma can be easily proved with basic computations. 
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Lemma 6.1. (Property of the convergence in measure) Let (f n )n, f C LP(£l) a sequence 
of functions such that 

• fn^-fin measure; 

• II /n 1 1 z,p (n) uniformly bounded. 
Then 

fn—tfin L q ({l), for every 1 < q < p. 

Another useful result in our proofs is a lemma concerning nonlinear monotone operators due 
to Brezis [BJ. 

Lemma 6.2. Let A be a maximal monotone operator on a Hilbert space H . Let Z n and W n 
be measurable functions from Q (a finite measure space) into H . Assume Z n — >■ Z a.e. on £1 
and W n -± W weakly in L 1 (Q;H). If W n (x) G A(Z n (x)) a.e. on O, then W(x) E A{Z(x)) 
a.e. on O. 

6.2 Regularity 

Concerning the regularity of the solution u of the DNLE, we refer for example to [13], [17] . 
[231. 



(Theorem 2.1 from [13j- inner Holder estimate) Let u be a weak solution of the DNLE. 
Then 

(6.1) u G C^ c p ' a {[0, T] x n) for some a G (0, 1). 
Moreover, for every cylinder Q' = [e,T] x Q' , £1' C £1, e > 0, we have 

(6.2) sup - < x, 
where a G (0, 1) and K > depend only on the T, f2' and the data. 

(Theorem 2.2. from |13j- Holder estimate up to the boundary) If 17 has regular 
boundary then 

u G C?J c p ' a ([0,T] x 0) for some a G (0, 1) 



and u satisfies an estimate similar to (6.2). 



6.3 Distance to the boundary function 

We collect some properties of the distance to the boundary function for which we refer to [TT] 
and [19]. Let d : 17 — > [0, +oo) be given by 

d(x) = dist(x, <9f2) = min{|x — z\ : z G 50}, 

where | • | is the Euclidean norm of M. d . In terms of d{x) we define the sets 

&i,r = {x £ £2 : d(x) > r}, 
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Q r = Q \ Qj t1 . = {i £ !] : d{x) < r}, 
and we remark that, for all small r > 0, 



dQi :r = {x £ Q : d(x) = r}. 

Lemma 6.3. (Properties of the distance to the boundary) Let ft C M> N be a bounded 
domain with boundary dQ of class C 2 . Then 

1. there is a constant £o G M + such that for every x G Q^ , there is a unique h(x) G dQ 
which realizes the distance 

d{x) = \x — h(x)\. 

Moreover, d(x) G C 2 (Q^ ), and for all r G [0,£o) the function H r : d(Q r ) H O, — > dQ, 
defined by H r {x) = h{x) is a homeomorphism. 

2. Function d(x) is Lipschitz with constant 1, i.e. 

\d(x) - d(y)\ <\x- y\. 

Moreover, 

(6.3) < c < |Vd(x)| < 1, Vx G 

and there exists a constant K > such that 

-K < dfjd(x) < K, \/x G , Vi, j = 1, N. 

Notice that this £o can be characterized as follows: 

£n = I min maxr : B r (x + rv) is tangent at in x}, 

xedn r>0 

where v is the inward unit normal at dQ in xq- We observe that d(xo + rv) = R and 
(6.4) a r C [J B r (y). 

y£dQ[ tT . 

Throughout the paper we have constantly used the notation £o with the properties stated 
above. 

7 Comments and open problems 

• In this paper we have discussed only the case slow diffusion case m(p — 1) > 1 and the 
quasilinear case m{p — 1) = 1. The fast diffusion case m(p — 1) < 1 produces different results 
and thus it needs different techniques. For this last case we mention the results of Savare and 
Vespri (|18j) about the asymptotic behaviour of the DNLE in the singular case. In that paper 
the authors prove the convergence to an asymptotic profile for a sequence of times t n — > T , T 
being the extinction time. The uniform convergence for all times and the rate of convergence 
in the fast diffusion case, for both DNLE and PLE, remain an open problem at this moment. 
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However, the fast diffusion regimes for the PME and the PLE have been much discussed in 
the literature, we mention [p] [TU]. 

• We presented only a formal description of the self similar solutions of the DNLE in the case 
m(p — 1) > 1. We do not offer a complete characterization of such functions since this beyond 
the purpose of our paper. The problem is interesting and it deserves a separate study itself. 

• Our result in the quasilinear case is not as sharp as the result in the degenerate case. Indeed, 
we only prove convergence in relative error. The problem of a rate of convergence similar to 
Theorem 1.1 is still open, except in the linear case m = 1, p = 2, where a representation as 
infinite series follows from the Fourier analysis of the solution. 

• For the Cauchy problem we mention [I] where the authors prove an L 1 -algebraic decay of 
the non-negative solution to a Barenblatt-type solution for the case j^jh < m < N^p-i) ' an< ^ 
they estimate its rate of convergence. 

• More general problems of this type can be considered by similar techniques. Let us mention 
the doubly nonlinear equation with mixed boundary conditions or p-Laplacian type equations 
with variable coefficients. 
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